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ABSTRACT 

^ ■ BP Psc is a remarkable emission-line field star that is orbited by a dusty disk 

and drives a parsec-scale system of jets. We report the detection by the Chandra 
X-ray Observatory of a weak X-ray point source coincident with the centroids 

P^ ■ of optical/IR and submillimeter continuum emission at BP Psc. As the star's 

photosphere is obscured throughout the visible and near- infrared, the Chandra X- 
ray source likely represents the first detection of BP Psc itself. The X-rays most 
likely originate with magnetic activity at BP Psc and hence can be attributed 
either to a stellar corona or to star-disk interactions. The log of the ratio of 
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Ph' X-ray to bolometric luminosity, log{Lx/Lhoi), lies in the range —5.8 to —4.2. 

O ■ This is smaller than log (Lx/L{,oi) ratios typical of low-mass, pre-main sequence 

'^ ■ stars, but is well within the log{Lx/Lhoi) range observed for rapidly-rotating 

^ ■ (FK Com-type) G giant stars. Hence, the Chandra results favor an exotic model 

wherein the disk/jet system of BP Psc is the result of its very recently engulfing 
a companion star or giant planet, as the primary star ascended the giant branch. 

^ ' Subject headings: circumstellar matter — stars: emission-line — stars: individual 

cn I (BP Psc) — X-rays: stars 

l> 

O 

o 

1. Introduction 

^ ■ Since its discovery by Stephenson (1986), the remarkable field star BP Psc (= StHa 

j^ ■ 202) has been classified as a classical T Tauri star (cTTS) on the basis of its strong Ha and 
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forbidden line emission. The discovery by Zuckerman et al. (2008; hereafter ZMS08) of an 
orbiting, dusty circumstellar disk and an enormous (pc-scale) system of highly collimated 
outflows (jets) would appear, at first glance, to provide strong support for such a cTTS 
classification. However, as noted by ZMS08, it is not at all clear that BP Psc is a pre-main 
sequence (pre-MS) star. 

First, although its jet system resembles those of cloud-embedded young stellar objects, 
BP Psc is located at high galactic latitude, far from any known star forming cloud. Fur- 
thermore, unlike the rare (though increasingly scrutinized) "isolated," cTTS systems within 
~ 100 pc (Torres et al. 2008) — the best-known examples being TW Hya (Kastner et al. 
1997; Zuckerman & Song 2004; and references therein) and V4046 Sgr (Kastner et al. 2008b; 
Rodriguez et al. 2010; and references therein) — BP Psc is not associated with any known 
sparse group of young stars (ZMS08). Second, its photospheric Li abundance appears to 
be anomalously weak for a K-type pre-MS star of age < 100 Myr, and its gravity-sensitive 
photospheric absorption lines also call into question such a classification (ZMS08). Third, 
the mm-wave molecular spectrum of its circumstellar disk — specifically, its weak CN and 
HCO"'" emission — is atypical of the disks of low- mass, pre-MS stars (Kastner et al. 2008a). 

In light of the first two considerations, ZMS08 proposed that BP Psc is most likely a first- 
ascent giant at a distance ~ 300 pc. They speculated that the BP Psc disk/jet system might 
therefore be the result of a recent, catastrophic interaction in which a low-mass companion 
star (or perhaps giant planet) had been consumed by BP Psc. Though somewhat exotic, a 
similar (companion engulfment) scenario may also apply to the first-ascent giant star TYC 
4144 329 2 (Melis et al. 2009): the post-MS nature of TYC 4144 329 2 is well-estabhshed, 
yet (like BP Psc) this star is evidently actively accreting gas from a dusty circumstellar disk. 

Low-mass, pre-MS stars are highly luminous X-ray sources, with typical (0.3-8.0 keV) X- 
ray luminosities log Lx [erg s~ ] ~ 29.5 and X-ray to bolometric luminosity ratios in the range 
\og{Lx/Lboi) ~ —4 to —3 (e.g., Preibisch & Feigelson 2005; Giidel et al. 2007a). Hence, X- 
ray observations are diagnostic of stellar youth. Furthermore, in rare cases, TTS jet systems 
may yield shocks energetic enough to produce X-rays, yielding insight into physical conditions 
within the jets (e.g., Kastner et al. 2005; Guedel et al. 2005, 2007b). Here, we report the 
results of Chandra X-ray Observatory observations of BP Psc intended to determine whether 
BP Psc and/or its jet system emit X-rays and, thereby, to better understand the nature of 
this unusual system. 
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2. Observations 

Chandra observed BP Psc in 2009 January, witli tlie Advanced CCD Imaging Spec- 
trometer (ACTS) as the detector. Due to spacecraft thermal constraints, the observation was 
obtained in two pieces of duration 10.15 and 65.4 ks (OBSIDs 8900 and 10856, respectively, 
on 2009 Jan. 12 and 13), for a total exposure time of 75.5 ks. Data were obtained with ACTS 
CCDs SI, S2, S3, S4, 12, and 13 active. Each CCD has a ~ 8' x 8' field of view with 0.49" 
pixels, and is sensitive to photons in the energy range 0.5-8.0 keV (with some additional 
sensitivity down to 0.3 keV, for back-illuminated CCD S3). The telescope was pointed such 
that BP Psc was positioned at the nominal aimpoint of CCD S3, ~ 2' from the S2-S3 chip 
boundary. The merged event data from the two exposures were subject to standard process- 
ing, filtering, and calibration (via Chandra X-ray Center (CXC) ACTS pipeline v7.6.1 and 
CALDB v3.5.1). 

We ran the CXC (CIAQ^ source detection tool celldetect on the filtered event list, 
restricted to CCDs S2 and S3 (the other active CCDs lie far off the best focal surface). We 
find 26 sources, including BP Psc itself (§3.1), were detected at a significance > 3cr (using 
Poisson statistics) on these two CCDs (Fig. [1], right). Four of these S2 and S3 sources lie 
within ~ 1" of DENIS near-IR sources; two of the four (BP Psc itself and an X-ray source 
detected on CCD S2; §3.3) are also 2MASS sources. Based on the ~ 10 X-ray sources that 
lie within ~ 1" of Guide Star Catalog stars (not including the BP Psc source), we estimate 
the absolute astrometric accuracy of the two-part Chandra pointing to be ~ 0.2". 



3. Results 

3.1. The X-ray source at BP Psc 

In Fig.[T]we display the wide-field Ha image originally published in ZMS08 alongside the 
broad-band (0.5-8.0 keV) Chandra/ ACIS image of the same field surrounding BP Psc. The 
comparison reveals that Chandrahas detected a faint, point-like X-ray source at the position 
of BP Psc. Source detection indicates a background-subtracted count rate of 0.24±0.07 ks~^ 
for this source, or ~ 18 source photons detected during the ~ 75 ks exposure. The BP Psc 
X-ray source lies at J2000 coordinates a =23:22:24.72, 6 = —02:13:41.66, with a positional 
uncertainty that is dominated by the estimated astrometric accuracy (±0.2"). Given the 
uncertainties, this X-ray position is precisely coincident with that of the 2MASS position of 
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BP Psc {a =23:22:24.70, 6 = -02:13:41.40) as well as with the centroid of submillimeter 
continuum emission from BP Psc (a =23:22:24.72±0.02, 6 = -02:13:41. 55±0.2; ZMS08). 
Analysis of photon arrival times (via a Kuiper test against a Poisson process; Press et al. 
1992) yields no evidence that the BP Psc source is variable, although we are unable to rule 
out source variability. 

Given the small number of photons detected, it is not reasonable to attempt to deduce 
the characteristic source temperature Tx and line-of-sight absorbing column Nh and, hence, 
intrinsic source flux Fx via standard X-ray spectral analysis. However, useful constraints 
on Nh and Fx can be obtained via the median energy of the detected source photons. 
Specifically, Feigelson et al. (2005) demonstrated that Orion Nebula Cluster (ONC) pre-MS 
star X-ray sources detected by Chandra/ ACIS-I display a relatively tight correlation between 
median X-ray energy and Nh- Getman et al. (2010) refined and further generalized the 
Feigelson et al. (2005) results, demonstrating (via spectral simulations) that the relationship 
between median X-ray energy and Nh should in principle apply to any source that can 
be modeled in terms of two-temperature optically thin thermal plasma emission suffering 
intervening absorption, and is applicable over a wide range of intrinsic source luminosity 
(i.e., for Lx from 10^'^ to 10^^ erg s— 1). As the Getman et al. analysis was carried out for 
sources detected on the front-illuminated ACTS CCDs (10-13), which have lower QE than 
CCD S3 at low energy (i.e., photon energies ~ 0.5 keV), application of the Getman et al. 
methods might result in an underestimate in the value of Nh obtained via the median photon 
energy of a source on S3. This is unlikely to be an important effect for a source as hard as 
that detected toward BP Psc in 2009, however, especially given the overall degradation in 
the soft (< 1 keV) X-ray sensitivity of ACIS in the decade since Chandra's launclo. 

To determine the median detected photon energy for the BP Psc source, we applied the 
ACIS Extract (AE) software package g (Broos et al. 2010) to extract and analyze the source 
photons. Doing so, we find the median background-subtracted photon energy of the BP Psc 
source in the 0.5-8.0 keV energy band is 3.47±0.33 keV. This median energy corresponds 
to an absorbing column log iV// [cm~^] = 22.9 ± 0.2, based on the "best model family" 
determined from the Getman et al. (2010) statistical spectral analysis method (we display 
the AE-extracted X-ray spectral energy distribution of BP Psc, overlaid with the Getman 
et al. "best model family" spectral model, in Fig. H]). The Getman et al. (2010) method 
yields an estimated apparent (measured) broad-band flux of Fx = (5.3 ± 1.4) x 10^^^ erg 



^http://cxc. harvard.edu/cal/memos/contam_memo.pdf 

^The ACIS Extract software package and User's Guide are available at 
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cm~^ s~^ (0.5-8.0 keV), corresponding to an intrinsic flux Fx = 8 x 10~^^ erg cm"^ s~^. The 
measured and intrinsic hard-band (2.0-8.0 keV) fluxes are (4.6 ± 1.4) x 10^^^ erg cm~^ s~^ 
and 1.4 x 10"^'^ erg cm"^ s~^, respectively. Although the statistical error on the broad-band 
intrinsic flux (-Fx) is small (~ 30%), the value of Fx is uncertain by almost a factor of 5 (and 
the hard-band intrinsic flux is uncertain by a factor ~2), due to systematic uncertainties in 
determining the generic spectral model that best describes the BP Psc source (see Getman 
et al. 2010, their Fig. 1). However, the foregoing estimates for apparent and intrinsic flux 
are consistent with those we obtained via simple simulations to match the BP Psc source 
ACIS-S3 count rate using absorbed thermal plasma models ranging in temperature from 
kTx = 1.0 to 3.0 keV — a range typical of coronal sources with median photon energies 
similar to that of BP Psc (Feigelson et al. 2005; also see §4) — assuming Nh = 10^^ cm~^. 

From the preceding results for Nh and Fx-, we estimate that the X-ray luminosity of the 
BP Psc source is Lx ~ 9 x 10^^(D/100)^ erg s~^ where D is the distance in pc. Adopting a 
bolometric luminosity Li,oi ~ 8.5 x 10^^(-D/100)^ erg s~^ for the star (ZMS08), and allowing 
for the (large) systematic errors in Fx obtained via the Getman et al. (2010) method, we 
conservatively conclude that log (Lx/Lhoi) lies in the range —5.8 to —4.2, with a best estimate 
onog{Lx/L,oi)^-5.0. 



3.2. Limits on diffuse X-ray emission 

As is apparent from inspection of Fig. 1, no diffuse X-ray emission is associated with 
the jets and bow-shock-like Herbig-Haro (HH) objects emanating from the position of BP 
Psc. This result is not particularly surprising, considering such emission should be very soft; 
indeed, there are few examples of shock-generated X-ray emission among pre-MS jet/HH 
systems, and the handful that have been detected to date generally display low X-ray lumi- 
nosities (e.g., HH 210 in the ONC, with observed Lx ~ 10^^ erg s~^; Grosso et al. 2006). 
We estimated an upper limit on the surface brightness of such diffuse X-ray emission from 
the BP Psc jet system by counting background events within regions centered on the Ha 
emission near the tips of the bow shocks seen in optical emission lines. The resulting limit 
is ~ 10~^^ erg cm^^ s^^ arcsec^^. This corresponds to an upper limit on emergent X-ray 
luminosity of Lx ~ 10^^(i5/100)^ erg s^^ for any compact (~ 1") X-ray-emitting knots in 
the BP Psc jet system. 



-6- 



3.3. Field X-ray sources 

If BP Psc is a nearby {D < 100 pc), "isolated" classical T Tauri star akin to TW 
Hya or V4046 Sgr (e.g., Kastner et al. 2008b and references therein), then it could be a 
member of a sparse group of young (age < 30 Myr) stars whose members are quite X-ray 
luminous relative to the field star population (Zuckerman & Song 2004; Torres et al. 2008; 
and references therein). If so, there would be a small chance that the ~ 0.04 square degree 
field imaged by Chandra (Fig. [1]) contains one or more low- mass members of this putative "BP 
Psc Association." Such young stars, being close to Earth, should be sufficiently luminous in 
X-rays as to be readily detectable in our ~ 75 ks Chandra exposure and should be bright 
near-lR sources (BP Psc itself, though occulted by its disk and not directly detected in the 
near-lR, has K ^ 7). 

As already noted, however, only one of the roughly dozen X-ray sources detected on 
CCD S3 — BP Psc itself — is also a 2MASS source. An additional dozen, relatively bright 
(> 50 count) X-ray sources are detected on S2 (the CCD abutting S3), although some of 
these sources suffer from poor image quality due to a combination of large off-axis angle 
and non-optimal detector geometry with respect to the focal surface of the Chandra mirror 
assembly Of these S2 sources only one, at a =23:22:41.4, (5=-02:21:23 (J2000) — the bright 
X-ray source at the extreme southeast corner of the image in the righthand panel of Fig. [T] — 
can be readily associated with a bright {V < 15) star (USNO-Bl 0876-0839011). Given its 
USNO and 2MASS catalog magnitudes {V = 14.09, R = 13.65, J = 12.36, H = 11.92, and 
K = 11.80) this X-ray-bright object is indeed a candidate young (age ~ 10 Myr), late- type 
(K3-4) star; but as such it would lie at a distance ~ 500 pc. Furthermore, its proper motion 
{Aa,AS = [48, 28] mas yr~^; USNO-Bl) is signicantly different from that measured for BP 
Psc ([44.4, -26.3] mas yr'^; Tycho-2). Hence USNO-Bl 0876-0839011 is not associated 
with BP Psc, whether or not BP Psc is a young star ai D < 100 pc. 



4. Discussion 

4.1. The X-ray source at BP Psc 

Our Chandra observation of BP Psc has yielded a detection of a faint X-ray point source 
that is coincident with the bright optical/near-lR and submillimeter emission associated with 
this unusual field star. The nature of this X-ray emission source is difficult to ascertain, given 
the small number of photons detected. However, the median energy (~ 3.5 keV) and, hence, 
the resulting inferred absorbing column density {Nh ~ 10^^ cm~^) appear typical of pre-MS 
star X-ray sources detected within dusty circumstellar disks that are viewed nearly edge-on 
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(Kastner et al. 2005). The large value we infer for Nh, which corresponds to a l^ band 
extinction of ~ 100 magnitudes (Vuong et al. 2003), is consistent with the fact that the 
photosphere of BP Psc is thus far only detected indirectly — via scattering off circumstellar 
dust — in the visible and near-IR (ZMS08). 

Given that the median detected photon energy is ~ 3.5 keV, it would appear that the 
X-rays detected from BP Psc are too hard to be attributed to internal shocks associated with 
its jet-launching regions — in contrast to cases in which soft (< 1 keV) X-rays have been 
detected at the bases of optical jets associated with certain classical T Tauri stars (Kastner 
et al. 2005; Guedel et al. 2005, 2007b). Indeed, the BP Psc X-ray source spectral energy 
distribution (Fig|2]) reveals no evidence for such a soft X-ray "excess." Instead, the X-ray 
photons we have detected would appear to be associated with magnetic reconnection activity 
near the star. Such magnetic activity, in turn, could be due either to a stellar corona or 
to star-disk interactions. In either case — given that the stellar photosphere is evidently 
obscured from direct view, even at near-IR wavelengths (ZMS08) — it would appear that 
the X-ray source at BP Psc likely represents the first direct detection of radiation from 
within a few stellar radii of BP Psc itself, and that the X-ray-based value of N^ we have 
obtained is therefore the only presently available measurement of the extinction toward the 
star {Ak ~ 10). 



4.2. Constraints on a pre-main vs. post-main sequence nature 

The values of X-ray luminosity and log (Lx/Lhoi) (§3.1) inferred from the Chandra 
detection would appear to place BP Psc at the extreme low end of the (well-determined) 
X-ray luminosity distributions of T Tauri stars in the ONC (Preibisch & Feigelson 2005) 
and in Taurus (Guedel et al. 2007a). Specifically, the value of Lx we estimate for BP Psc, 
Lx ~ 10^^ erg s~\ places BP Psc well below the median value of Lx for 0.4 — 2.OM0 ONC 
stars in the 3-20 Myr age range (the relevant range if BP Psc is pre-main sequence; ZMS08). 
Only - 2% of ONC stars display Lx < 3 x lO^^ erg s"^ (Preibisch & Feigelson 2005). 

In contrast, the inferred range of log (Lx/Lhoi), —5.8 to —4.2, would place BP Psc well 
within the (rather broad) \og{Lx/Lboi) distribution of the class of rapidly rotating (FK 
Com-like) G giants. The dozen or so members of this class range from log (Lx/Lhoi) = —6.4 
to —3.1, with a median of \og{Lx/Lboi) = —4.7 (Gondoin 2005). Hence, the Chandra 
identification of a weak X-ray source at the position of BP Psc appears to throw some 
weight behind a post-main sequence classification for this star — and supports a model 
in which the detected X-rays most likely arise from coronal activity, as a consequence of 
its rapid rotation rate (see ZMS08). The large angular momentum required to sustain such 



rapid rotation, in turn, would seem to be consistent with the notion that BP Psc has recently 
consumed a companion star (as proposed by ZMS08), and is therefore an FK Com-type star 
"in the making" (as speculated by Kastner et al. 2008a). 



5. Summary and Conclusions 

We have used the Chandra X-ray Observatory to detect a weak X-ray point source that 
is coincident with the centroids of optical/IR and submillimeter continuum emission from 
the disk-enshrouded, jet-driving field star BP Psc. As the star's photosphere is obscured 
throughout the visible and near-infrared, the Chandra X-ray source appears to represent the 
first detection of BP Psc itself. The Chandra X-ray image yields no detection of the BP Psc 
jet system, nor evidence that BP Psc is associated with a young stellar group or cluster. 

From the median energy of the detected BP Psc source photons (~ 3.5 keV), we infer 
an absorbing column log Nh [cm~^] = 22.9 ± 0.2 and an X-ray luminosity Lx ~ 10^^ erg 
s^^), with \og{Lx/Lboi) in the range —5.8 to —4.2. The rather large median X-ray energy 
indicates that the X-rays most likely originate with magnetic activity, as opposed to shocks 
near the star, and hence are produced either in a stellar corona or in star-disk magnetospheric 
interaction regions. The inferred values of Lx and Lx/Lhoi are 1-2 orders of magnitude 
smaller than those typical of low-mass, pre-main sequence stars, but are well within the 
range observed for rapidly-rotating (FK Com-type) G giant stars. Hence, the Chandra 
results favor an exotic model wherein the disk/jet system of BP Psc is the result of its very 
recently engulfing a companion star or giant planet, as the primary star ascended the giant 
branch. 
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Fig. 1. — Ze/t; wide-field Ha image of BP Psc and surrounding field, showing the jet system 
associated with BP Psc (adapted from Zuckerman et al. 2008). Right: broad-band (0.5- 
8.0 keV) 75.5 ks exposure Chandra/ ACIS image of the same field. Only X-rays detected 
by ACIS CCDs S2 (lower region) and S3 (upper region) are displayed; the field of view is 
8' X 16'. The circles on each image indicate the positions of Chandra X-ray sources, with 
radius proportional to signal-to-noise ratio (specifically, radii [arcsec] are set to 2x the source 
signal-to- noise ratio; the smallest circles correspond to 3o" detections). In each image, north 
is up and east is to the left, and the position of BP Psc is indicated by an arrow. 
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Fig. 2. — Spectral energy distribution of the X-ray source at BP Psc (crosses) overlaid with 
the "best model family" spectral model (histogram; see Getman et al. 2010). The parameters 
describing this (two-component thermal plasma) model are kTi = 0.8 keV and kT2 = 1.8 
keV, with emission measure ratio EMi/EM2= 1.0, metal abundance Z = 0.3 (relative to 
solar), and intervening absorbing column log Nh [cm~^] = 22.9. 



